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Abstract 


It is essential to know the space-time structure of the nonlocal vacuum 
condensates for application to medium energy processes. Using the Dyson- 
Schwinger formalism in the rainbow approximation for the quark propagator, 
we study the nonlocal quark condensate and model forms for the nonper- 
turbative gluon propagator constrained by fits to local condensates and deep 
inelastic scattering with nucleon targets. 
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I. INTRODUCTION 


It has long been known that chiral symmetry breaking reqnires a nonpertnrbative qnark 
propagator with non vanishing vacnnm matrix elements of normal ordered products of quark 
helds m, called quark condensates. These condensates would vanish in a perturbative 
vacuum but do not vanish in the QCD vacuum, and are of central importance both for the 
structure of hadronic matter and for the study of the early universe chiral phase transition. 
A systematic treatment of hadronic masses can be carried out |^|-Q using operator product 
expansions in terms of vacuum matrix elements of local operators, the vacuum condensates, 
whose phenomenological values have been conhrmed in lattice gauge calculations |^. As 
was discussed in early work on the magnetic dipole moments of nucleons 0 , for application 
to form factors and transition matrix elements in the low to medium momentum transfer 
region the operator product expansion cannot be used, since long distance properties of 
nonlocal operators must be treated. 

One approach to this problem of treating bilocal operators has been the use of nonlocal 
condensates, which have been introduced to represent the bilocal vacuum matrix elements 
needed for the pion wave function and pion form factor for low to medium momentum 
transfer. In this method one does not carry out an O.P.E. for the vacuum matrix elements of 
the bilocal operators, but introduces new phenomenological functions needed to characterize 
the space-time structure of the nonlocal condensates. Both the forms of these functions and 
the parameters are found by hts to experiment as well as considerations of analyticity. E.g., 
in a study of parton distribution functions 0 the space-time scale of a nonlocal condensate 
was determined by a £t of a monopole form in space-time to experimental data. On the 
other hand, in a recent use of nonlocal condensates to determine the values of vacuum 
susceptibilities |T^ , which characterize the nonpertnrbative quark propagation in an external 
held o, it was found that the monopole form did not have suitable analytic properties, 
and a space-time dipole form was used to £t the low-x parton data. Although a satisfactory 
£t to the phenomenological pion susceptibility 0 was found, it is a good example of 
the importance of determining the structure of the nonlocal condensates for application to 
transition matrix elements over a wide range of momentum transfer. 

It is the goal of the present work to study the form of the nonlocal quark condensate using 
the the QCD Dyson-Schwinger [D-S] equations |T3| , p!^ . Using the bare gluon-quark vertex, 
dehned as the rainbow approximation, the nonpertnrbative (dressed) quark propagator is 
determined self consistently with a model for the nonpertnrbative (dressed) gluon propagator 
D“^(g). A comprehensive review of this type of model is given in Ref. []^ It has been 
shown that in the rainbow approximation the value of the quark condensate |]l6ipT| and the 
mixed quark condensate can be obtained with suitable choices of the gluon propagator, 
which also provides constraints for the present work. The gluon condensate within this 

m 


approach has been studied in Ref. 
approximat ions 


Other studies use the D-S formalism with different 


201 to attempt to determine the nonpertnrbative quark condensates. 


For hadronic properties such as the elastic and transition form factors one needs the 
information equivalent to the bound-state Bethe-Salpeter [B-S] equation. It has also been 
shown that the rainbow D-S model is consistent with nontopological chiral quark models 
g and low-momentum transfer meson form factors m- J^or a treatment of form factors 


over an extended range of momentum transfer light-cone B-S studies of the pion form factor 
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have shown that nonperturbative as well as perturbative two-quark propagators are 
needed, and that in a QCD treatment four-quark matrix elements are also required. We 
plan a study of nonlocal nonperturbative four-quark matrix elements within this approach 
in the near future. 

In the present work we use the rainbow Dyson-Schwinger equation to investigate the 
forms of the nonlocal quark condensate as well as the gluon propagators. Using the same 
phenomenological gluon propagators as were used in previous studies of the local condensates 
we hnd that the dipole form with the parameter close to the one found from hts to 

,lTni 


the sea-quark distribution 


can be obtained. 


II. NONLOCAL QUARK CONDENSATE FROM AN EFFECTIVE 
INTERACTION IN THE DYSON-SCHWINGER APPROACH 


The quark propagator is dehned by 

5,(a;) = <0|T[g(a;)g(0)]|0>, 


( 1 ) 


where q{x) is the quark held and T the time-ordering operator. For the physical vacuum 
the quark propagator Sq{x) has a perturbative and a nonperturbative part. In the case of 
vanishing current quark masses (mo = 0) one can write 

Sqix) = S^^ix) + S^^ix), 


(x) = 


1 7 • X 


27r2 x'^ 


S. 


NP/ 


X = 


12 


(<: q(x)q(0) :> +x^ <: q(x)^^q(0) :>) 


( 2 ) 


It should be stressed that normal-ordered products, and therefore do not vanish in the 
nonperturbative vacuum. For short distances, the O.P.E. for the scalar part of S^^{x) gives 


X 


<: q{x)q{0) :>=<: g(0)g(0) :> -^ < 0| : g(0)a ■ G(0)g(0) : |0 > + ..., (3) 

in which the local operators of the expansion are the quark condensate, the mixed conden¬ 
sate, and so forth. 


In Ref. it is shown that with a choice of nonlocal condensate 
< 0| : q{x)q{0) : |0 >= g(x^) < 0| : g(0)g(0) : |0 >, 


( 4 ) 


with 


g{x^) = 


f(a) = 


(1 + 


10 


daf{a)e 






( 5 ) 
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and = 0.15... 0.20GeV^, one can fit the low-x quark distributions and also the pion 
susceptibility. 

In the Dyson-Schwinger formalism is related to the quark self-energy, S, by 

Sq{p)~^ = il ■P + T.ip). ( 6 ) 

Using the bare quark gluon vertex (the rainbow approximation), r|)(g,p) = 7i/^, S(p) 
satishes the rainbow Dyson-Schwinger equation ||14|| : 

s(p) = / - 9)7^ fs5(9)7.y (7) 

with the gluon propagator, A“ the color SU{3) matrix. In Euclidean space one can 

write 


Sqip) =i-fpA{p) + B{p) 


( 8 ) 


The choice of the Landau gauge for the gluon propagator 


Dtiq) = r 




fiiy 


D{q^ 


(9) 


leads to the set of coupled integral equations 




B{p^) = ^gl J 


A{q^ 


q2j^2(^q2) _|_ 

B{q^) 

q2j^2(^q2) _|_ B{q^Y 


p ■ q + 2 


{p-q-q^){p^ -p-q) 


{p-qf 


( 10 ) 


Mesonic bound states can be studied within this framework by solving the ladder Bethe 
Salpeter equations for the corresponding qq bound states. Various mesonic properties have 
been studied in Refs. [0. In Ref. a detailed investigation of the low energy sector was 
performed by deriving the general form of the effective chiral action for the SU(3) Goldstone 
bosons and determining and most of the chiral low energy coefficients Li, which, in turn, 
determines the physics of the tt, K and p mesons at low energies . 

Because the form of the gluon propagator D{s) in the IR region is unknown, we must 
use model forms as input. Our model ansatz is 


gs^D{s) = a , 


( 11 ) 


which determines the quark-quark interaction through a strength parameter y and a range 
parameter A. Its form is inspired by the 5 function ansatz of Ref. [M|, which it approaches 
for A —> 0. 

The nonlocal quark condensate <: ^(x)g(O) :> is then given by the scalar part of the 
Fourier transformed inverse quark propagator: 
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<: q{x)q{0) :> = (-4A^c) 


12 


d^p 


B{p^ 


{27r)^ p‘^A‘^{p‘^) + B‘^(jP') 


^ipx 


dss 


B{s) 


Ji 


Jo sA2(s) + 

At X = 0 the expression for the local condensate <: qq :> is recovered: 


<; qq ;>= (-) 


12 


^OO 


dss 


B{s) 


16^2 Jo s/12(s) + B^s) 


( 12 ) 


( 13 ) 


The nonlocality gi^x"^) can be obtained immediately by dividing ( 0 ) through ( 0 )- 

Because the quark-quark interaction dehned by (^) has a hnite range in momentum 
space the momentum integrals in (0 and ( 0 ) are hnite. Our analysis ignores effects from 
hard gluonic radiative corrections to the condensates which are connected to a possible 
change of the renormalization scale p at which the condensates are dehned. Those ehects 
are of minor importance for our study of nonperturbative ehects in the low and medium en¬ 
ergy regions. It should be stressed in this context that our interaction is not renormalizable 
because we are using the bare quark gluon vertex. Therefore, instead of our condensates de¬ 
pending logarithmically on the renormalization scale p [^], the scale at which a condensate 


is dehned in our approach is a typical hadronic scale, which is implicitly determined by the 
model gluon propagator g^D{s) and the solutions of the D-S equations (0). The situation 
is very similar to the determination of vacuum condensates in the instanton liquid model 
where the scale is set by the inverse instanton size. 

In order to check the sensitivity of our results on the model gluon 2 point function (pH]) 
we try various sets of parameters y and A and investigate the dependence of the function 
for these forms. We solve the set of integral equations (pTf ) self consistently for a given 
model form for gJ^D^q^) obtaining the quark propagator functions and B{p‘^), which, 

in turn will allow us to calculate g{x^) from (|T^). The result can then be compared to the 
dipole fit of Ref. lOT with = 0.15... 0.20GeV^. The parameter sets we are using are 


Set 1: A = 2.0 * lO’^GeV^ ; y = l.dOGeV 

Set 2 : A = 1.0 * lO-^GeV^ ; y = 1.56GeV 

Set 3: A = 2.0 * lO’^GeV^ ; y = 1.58GeV. 


(14) 


These parameters have been chosen so that they reproduce the correct value for the pion 
decay constant in the chiral limit = 88 MeV. Moreover the values of the chiral low 
energy coefficients Lj are compatible with the phenomenological values in both cases. 
Following Ref. one hnds : Li = 0.86, L 3 = —4.53, L 5 = 0.78, Lg = 0.84 (*10“^) for 
Set 1; Li = 0.84, Lg = -4.48, L 5 = 0.88, Lg = 0.84 (*10-3) for Set 2 and Li = 0.83, 
Lg = -4.42, Ls = 0.92, Lg = 0.78 (*10-3) for Set 3. 


Many of the works in Refs. 


have used a model ansatz for the gluon propagator 


DtM) = 


ab ) 


( 16 ) 


which is often referred to as the Feynman-like gauge. It is however not identical to the 
Feynman gauge QGD, in which the dressed gluon propagator would have different longitu¬ 
dinal and transverse components. Therefore the ansatz (^) should be regarded merely as a 
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model form for the gluon 2 point function. For our purpose it is, however, interesting to ask 
if and how the nonlocal quark condensate depends on choosing either Landau gauge (|^) or 
the Feynman-like gauge (pISl). Therefore we perform the calculation for another parameter 
set: 


Set 4: A = 2.0 * IQ-^GeV^ ; x = 1.23GeV, 


(16) 


while using ([T5|) instead of for the gluon propagator. Set 4 has the same range parameter 
A than Set 1. The strength parameter y is slightly smaller in order to obtain the correct 
value of Jtt = 88 MeV. The chiral low energy coefficients are: Li = 0.85, L 3 = —4.46, 
L 5 = 0.82, Lg = 0.93 (*10“^), values rather close to those obtained with Set 1. 


III. RESULTS AND DISCUSSION 


Fig.^ shows the results for g{x‘^) for the four parameter sets and compares with the dipole 
£t of Ref. with = 0.20GeV^ (solid line). The result of Ref. is best reached for 
a gluon propagator with a small range parameter A = 0.002GeV^ in the infrared. Larger 
values for the width parameter A lead to stronger deviations from the form of Ref. |T0[ . 

By comparing the curves for Set 1 and Set 4 we can demonstrate that the nonlocal 
condensate is very robust with respect to using Landau gauge (H) or the Feynman-like 
gauge (|^). The change between the two forms for the gluon 2 point function can be easily 
made up by a slight readjustment of the parameters of the IR model ansatz without any 
signihcant change of the hnal result. 

We conclude that the D-S formalism is a valuable tool for the study of the nonlocal quark 
condensate, and expect that the B-S formalism will also prove to be useful for the study of 
the nonlocal four-quark condensates, which provide nonperturbative QGD effects for hadron 
couplings and form factors. 
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FIGURES 


FIG. 1. The non local quark condensate g{x) =<: q{x)q{0) :> / <: q{0)q{0) :> for the four 
sets of model gluon propagators mentioned in the text compared with the dipole fit of Ref. |1C]. 
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